Journal  of  Power  Sou 


;  249  (2014)  219-230 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


Diffusive  relaxation  of  Li  in  particles  of  silicon  oxycarbide  measured 
by  galvanostatic  titrations 

P.D.  Weidman,  Dongjoon  Ahn,  R.  Raj 

Department  of  Mechanical  Engineering,  University  of  Colorado  at  Boulder,  Boulder,  CO  80309-0427,  USA 


& 


CrossMark 


HIGHLIGHTS 


•  We  show  that  diffusive  relaxation  during  galvanostatic  titrations  embodies  one  principal  time  constant. 

•  The  time  constant  is  independent  of  the  rate  of  injection  of  Li  just  prior  to  the  open  circuit  galvanostatic  relaxation. 

•  The  time  constant  is  related  explicitly  to  the  particle  size  and  to  the  chemical  diffusivity  of  Li  in  the  active  material. 


ART 


C  L  E  INFO 


A  B  S  T  R 


C  T 


Article  history: 

Received  14  June  2013 
Received  in  revised  form 
18  September  2013 
Accepted  21  September  2013 
Available  online  2  October  2013 


Keywords: 

Lithium  ion 
Batteries 

Chemical  diffusion 
GITT 

Silicon-oxycarbide 


Experiments  with  silicon-oxycarbide,  an  anode  material  for  Li-ion  batteries,  are  compared  with  a  rigorous 
solution  for  time-dependent  diffusion  profiles  in  spherical  particles,  which  is  supplemented  by  finite 
element  analyses  of  cuboid  shaped  particles.  In  this  way  a  value  for  the  chemical  diffusivity  D  1.8  x  10  18 
-4.2  x  10  18  m2  s  1  across  the  entire  range  of  the  state-of-charge  is  obtained.  The  method  consists  of 
galvanostatic  titrations,  where  a  constant  current  injection  is  followed  by  the  measurement  of  the  galvanic 
potential  in  open  circuit  relaxation,  at  various  values  of  the  state-of-charge.  Comparison  with  theory  shows 
that  the  relaxation  time  varies  with  the  state-of-charge  because  the  particles  have  variable  size:  small 
particles  have  short  relaxations  and  dominate  the  low  values  of  the  state-of-charge.  The  theoretical  results 
have  additional  notable  features.  We  show  that  the  solution  for  the  spherical  particles  can  be  represented  to 
a  very  good  approximation  by  a  single  time  constant  of  exponential  relaxation.  Its  value  is  significantly 
slower  than  the  results  for  one-dimensional  solutions.  The  relaxation  time  is  independent  of  the  injection 
current  and  the  state-of-charge,  imparting  significant  generality  to  the  results. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  diffusivity  of  lithium  into  anode  and  cathode  materials 
controls  the  rate  at  which  batteries  can  be  charged  and  discharged. 
It  is  therefore  a  fundamental  parameter  for  materials-selection  in 
battery  design. 

The  effective  rate  of  Li  insertion  or  extraction  is  related  to  the 
physical  length  scale  of  the  electrode  material,  A,  and  to  the  coef¬ 
ficient  of  chemical  diffusion  of  Li  in  the  host  material,  which  we 
write  as  D.  From  dimensional  arguments  the  typical  time  for 
diffusion,  t,  across  this  length  scale,  can  be  normalized  according  to 
the  scaling  law 


A2  =  mDt  (1) 
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where  m  is  a  dimensionless  parameter  whose  value  can  vary 
significantly  depending  on  the  geometry  of  the  diffusion  field.  A 
key  objective  of  the  present  work  is  to  provide  exact  descriptions 
for  the  diffusion  of  Li  into  spherical  particles  in  ways  that  can  be 
easily  compared  with  experiment  in  order  to  elicit  values  for  the 
chemical  diffusivity  of  Li.  The  experiments  are  carried  out  on 
silicon-oxycarbide,  an  anode  material  for  Li-ion  batteries. 

The  architecture  of  a  typical  anode  is  illustrated  in  Fig.  1  [1  ].  The 
anode  particles  are  connected  electrically  to  the  current  collector 
via  a  conducting  agent  (typically  carbon  black).  The  voltage  is 
measured  across  the  half-cell  with  pure  Li  metal  as  the  reference 
electrode.  In  the  open-circuit  condition  when  no  current  flows 
through  the  cell,  this  voltage  is  related  to  the  chemical  activity  of  Li 
on  the  surface  of  the  particles,  aLi,  by  the  Nernst  equation 

V  =  *rin(au).  (2) 


0378-7753/$  —  see  front  matter  ©  2013  Elsevier  B.V.  All  rights  reserved. 
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where  F  is  the  Faraday  constant.  Since  activity  of  lithium  in  the  host 
must  be  less  than  that  of  pure  metal,  au  <  1,  and  therefore  V  <  0. 
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Nomenclature 

a,  r  the  radius  of  the  electrode  particle,  and  the  radial 
coordinate  from  the  center  of  the  particle. 
c(r,t)  the  spatial  and  time-dependent  concentration 

(moles  per  unit  volume)  of  the  diffusing  species. 
c(t)  the  volume-averaged  charge  state,  i.e.  the  total 

charge  (in  moles)  injected  into  the  particle  divided 
by  the  volume  of  the  particle. 
ca(t),  c„(t)  concentration  at  the  surface  of  the  particle. 
t  real  time 

t  dimensionless  time 

Jo  injection  rate  written  as  flux  per  unit  area  entering 

the  particle 

DF  Fickian  diffusion  coefficient 

D  chemical  diffusion  coefficient 

to.  to  period  of  current  injection  in  real  and 
dimensionless  time 

t*.  t*  real  and  dimensionless  time  measured  from  the 
start  of  open  circuit  relaxation 
At*  time  constant  for  the  exponential  relaxation  of 

surface  concentration  during  open  circuit 
Qu  Li  capacity  in  mAh  g_1 

SOC  the  state-of-charge  defined  as  the  fraction  of  the  full 

capacity  of  the  material  achieved  at  a  given  current 
density 

a  a  nondimensional  parameter  for  calculating  the 

relationship  between  Li  concentration  and  open 
circuit  voltage. 

av  a  modified  version  of  a  in  units  of  V. 


(Current  Flow  Direction) 


The  relationship  between  concentration  and  the  activity  of  Li  is 
discussed  in  Section  4. 

The  diffusivity  of  Li  in  electrode  materials  is  often  measured  in 
one  of  two  ways,  known  as  PITT  [2]  and  GITT  [3],  In  PITT  a  step 
voltage  is  applied,  and  the  current  transient  is  measured.  In  GITT  a 
short  pulse  of  current  is  applied,  followed  by  the  measurement  of 
the  voltage  response  in  open-circuit.  The  relaxation  time  is 
substituted  into  an  equation  having  the  form  given  in  Eq.  ( 1 )  for  an 
estimate  of  the  diffusivity.  The  GITT  method  is  expected  to  be  more 
reliable  for  the  determination  of  the  relaxation  since  it  avoids  errors 
stemming  from  interfacial  polarizations  and  resistive  potential 
drops  across  the  electrolyte  [4,5], 

The  theory  and  experiment  presented  here  are  concerned  with 
galvanostatic  titrations  for  particles  of  spherical  geometry  and 
cuboid  shapes.  The  parameters  are  normalized  in  order  to  obtain 
self-consistent  values  for  the  coefficient  of  chemical  diffusion  from 
experiments  carried  out  for  different  values  of  injection  currents, 
states-of-charge,  and  particle  size. 

The  following  section  begins  with  the  theoretical  results,  fol¬ 
lowed  by  experimental  measurements  which  are  then  compared 
with  theory. 


2.  Analysis 


The  diffusion  transport  of  Li  into  electrode  materials  is  often 
analyzed  with  Fick’s  first  law,  where  the  flux  is  equal  to  the  product 
of  the  diffusion  coefficient  and  the  concentration  gradient  [6],  We 
call  this  the  Fickian  diffusion  coefficient,  DF.  However,  in  situations 
that  involve  chemical  reactions  it  is  more  appropriate  to  consider 
the  gradient  of  the  chemical  potential  as  the  driving  force  for 
diffusion;  in  this  instance  the  chemical  diffusion  coefficient,  D,  is 
the  measure  of  diffusivity.  The  two  diffusion  coefficients  are  related 
by  the  change  in  the  chemical  activity  with  concentration.  In  the 
sequel  open-circuit-relaxation  measurements  are  used  to  deter¬ 
mine  the  value  for  this  factor. 

In  the  next  section  two  specific  diffusion  problems  are  solved. 
First  the  species  (Li)  is  inserted  into  the  particle  at  a  constant  cur¬ 
rent,  and  the  resulting  time-dependent  concentration  profile  of  Li 
in  the  particle  is  determined.  In  the  next  step,  the  current  is  stopped 
(the  open-circuit  condition)  and  the  relaxation  of  the  voltage  pro¬ 
file  with  time  is  analyzed.  Subsequently,  these  results  will  be  used 
as  the  basis  for  interpreting  experiments  of  Li  intercalation  into 
particles  of  silicon-oxycarbide. 


2.1.  Description  of  the  analysis 

The  problems  solved  below  are  described  in  Fig.  2.  Under  con¬ 
stant  current,  as  shown  on  the  left,  the  surface  concentration  rises 
monotonically  with  time.  Under  open  circuit,  shown  on  the  right,  the 
surface  concentration  decreases  with  time  as  the  concentration  re¬ 
laxes,  eventually  reaching  a  uniform  value  throughout  the  particle. 


2.1.1.  The  governing  equation 

The  governing  equation  for  time  varying  diffusion  that  enforces 
mass  balance  and  continuity  of  concentration  in  space  and  time,  for 
a  spherical  geometry  is  given  by 


dc(r,  t)  fd2c  2dc\ 

^^  =  DF(var2  +  7arJ- 


(3) 


2.1.2.  Normalization  of  variables 

The  volume-averaged  concentration  in  the  particle  is  readily 
obtained  as  the  product  of  the  injection  rate  into  the  particle  (equal 
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Fig.  2.  Schematics  of  concentration  profiles  of  Li  in  a  particle  during  insertion  at  a 
constant  current,  followed  by  diffusional  smoothing  in  open  circuit. 


to  the  flux  times  the  particle  surface  area)  multiplied  by  the 
duration  of  injection  divided  by  the  volume  of  the  particle.  These 
substitutions  lead  to 

m  »  (4) 


The  solution  to  the  governing  equation  subject  to  the  boundary 
conditions  specified  in  Eq.  (8)  is  given  in  Appendix  A.  We  thereby 
obtain  the  concentration  at  the  surface  of  the  particle,  expressed  in 
dimensionless  time  as 


c“w--Br'+w(jr24-ir| 


where  the  eigenvalues  fa  satisfy 


O) 


fa  =  tan  fa.  (10) 

In  order  to  understand  the  behavior  predicted  by  Eq.  (9),  it  is 
useful  to  write  it  in  the  dimensionless  form 


Cq(t) 

(*) 


The  asymptotic  behavior  of  Eq.  (11)  as  x  — ► 


(ID 


is  given  by 


+  X  as  t-°°.  (12) 


The  surface  concentration  increases  linearly  with  time  at  large 
time,  with  1/15  intercept  on  the  ordinate. 

A  reasonable  approximation  to  Eq.  ( 11 )  is  obtained  by  expanding 
the  sum  and  considering  only  the  first  term.  The  first  three  eigen¬ 
values  are  fa  =  4.5,  fa  =  7.7,  and  fa  =  10.9.  Since  terms  in  the 
expansion  decay  exponentially  with  fa,  only  the  first  term  would 
make  a  significant  contribution  leading  to  the  result 


Cfl(t) 

0*) 


The  analysis  leads,  naturally,  to  normalization  of  time,  t,  given 
by 


_  Dpt 

The  dimensionless  radial  coordinate  is 


(5) 


R 


a ' 


(6) 


Combining  Eqs.  (4)— (6)  gives  the  average  concentration  in 
terms  of  normalized  time  as 


c(t)  = 


3q/o 

DF  • 


(7) 


Note  that  we  analyze  the  problem  in  the  real-time  domain.  An 
alternative  approach  with  Laplace  transforms  has  been  used  by 
White  and  co-workers  [6], 


2.2.  Results  for  constant  current  injection 

In  this  problem  we  seek  the  change  in  the  surface  concentration 
with  time  when  a  constant  current  is  injected  into  the  particle 
which  implies  the  boundary  and  initial  conditions 

Jo  -  X3F9C|^|r_a,  =  0-  c(r>°)  =  0  (8) 

The  first  condition  specifies  that  the  radial  gradient  of  the 
concentration  at  the  surface  is  equal  to  the  injection  rate.  The 
second  condition  arises  from  the  spherical  symmetry  of  the  prob¬ 
lem,  and  the  last  relation  ensures  that  the  initial  concentration  is 
zero  everywhere  within  the  particle. 


+  B(,-k20')-  <13> 

The  exact  result  from  Eq.  (11),  with  the  asymptotic  override,  is 
shown  in  Fig.  3.  The  plots  show  two  regimes  of  behavior:  the 
transient  regime,  which  has  duration  of  t  =  0.05,  followed  by  a 
steady  state  regime.  The  transient,  which  we  call  At;,  with  the 
subscript  standing  for  the  constant  current  case,  is  related  to  the 
diffusion  coefficient  via  Eq.  (5) 


Ca(T) 

m 


0.05  a2 


(14) 


Plots  for  Eq.  ( 11 )  that  include  only  the  first  term,  fa,  two  terms  fa, 
fa,  and  three  terms  fa,  fa,  fa,  are  shown  in  Fig.  3.  The  result  from 


T 


Fig.  3.  Plots  of  Eq.  ( 1 1 )  showing  the  change  in  surface  concentration  with  time.  The  full 
plot,  as  well  as  the  asymptote  and  plots  for  one,  two  and  three  terms  from  the  sum¬ 
mation  of  the  exponential  terms  are  shown. 
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Eq.  (13),  which  includes  only  the  first  term,  gives  a  very  good 
approximation  to  the  exact  result  for  all  x  >  0.03. 

In  the  second  problem  considered  below,  the  current  is  stopped 
and  the  diffusing  species  are  allowed  to  relax  towards  a  uniform 
concentration  within  the  particle. 


2.3.  Analysis  for  open  circuit  relaxation 


In  this  section  we  analyze  the  problem  where  the  current  is 
stopped  at  t  =  to,  and  the  concentration  within  the  particle  is 
allowed  to  relax  by  diffusion,  eventually  leading  to  a  uniform 
concentration.  This  uniform  concentration  is  given  by  Eq.  (4), 
evaluated  at  time  to. 

To  track  the  time  in  the  open  circuit  state,  we  introduce  the  new 
time  variable,  t*,  defined  as 

t*  =  t-t0  (15) 

where  t*  begins  after  the  current  is  stopped.  The  normalized  time 
variables,  x0  and  r*  are  defined  in  the  same  way  as  is  x  in  Eq.  (5). 

This  problem  is  solved  in  Appendix  B,  with  the  boundary 
condition 


=0 

9r  lr=0 

and  with  the  t*  =  0  initial  condition 

m  3 j/ptb  3J0a  Jor2  2J0a  ^  sin Qg„JQ 

a  10Df  2 aDF  DF  ^  /j'^Rsin  /?„ 


(16) 


(17) 


where 


(18) 


The  analysis  gives  the  following  evolution  of  surface  concen¬ 
tration  with  time 


* 


2  ^ 
=  t°+3  Y, 

J  l 


(19) 


This  can  be  simplified  to  a  very  good  approximation  by 
considering  only  the  first  term  in  the  expansion,  viz. 


The  result  in  Eq.  (21)  is  particularly  relevant  to  experimental 
measurement.  It  predicts  that  the  surface  concentration  will  decay 
exponentially  at  the  same  rate  regardless  of  the  magnitude  and  the 
duration  of  the  current  injection. 

In  the  following  section  the  above  results  are  applied  to  elec¬ 
trochemical  data  from  half-cells  where  anodes  are  made  from 
particles  of  silicon-oxycarbide. 

3.  Experiments 

3.1.  The  method 

The  experiments  were  carried  out  with  electrochemical  inser¬ 
tion  of  Li  into  anodes  made  from  particles  of  amorphous  silicon- 
oxycarbide,  employing  Li  metal  as  the  counter  electrode.  These 
materials  are  made  by  controlled  pyrolysis  of  siloxanes,  which  are 
silicon  based  polymers.  They  are  amorphous.  Their  structure  con¬ 
sists  of  a  stochastic  network  of  graphene,  with  an  average  domain 
size  of  ~2  nm.  The  intradomain  spaces  are  filled  with  SiC>4  tetra- 
hedra,  but  the  interfaces  with  graphene  also  contain  mixed  bond 
tetrahedra  constructed  from  Si— C— 0  [7], 

The  procedure  for  the  preparation  of  the  powders  and  the  as¬ 
sembly  of  the  half  cells  are  described  in  Ref.  [  1  ].  The  composition  of 
the  amorphous  material  for  the  present  study  was  SiCi.gsOo.ss-  The 
average  size  of  the  particles  in  the  powder  was  ~  5  pm.  The  elec¬ 
trode  was  prepared  by  mixing  SiCO  powders  with  carbon  black  and 
a  binder  into  slurry,  which  was  screen  printed  on  a  copper  foil  with 
a  doctor  blade,  followed  by  drying  and  pressing.  The  thickness 
of  the  electrode  material  on  copper  produced  in  this  way  was  50— 
70  pm.  Lithium  foil  was  used  as  the  counter  and  reference  elec¬ 
trode.  The  electrodes  were  sandwiched  with  microporous  poly¬ 
propylene  (Celgard,  USA)  and  filled  with  the  electrolyte  consisting 
of  a  solution  of  1  M  LiPF6  and  diethyl  carbonate  (volume  ratio  1:1, 
Ferro  Corporation,  USA).  The  coin  cells  were  assembled,  crimped 
and  closed  in  an  Argon  filled  glove  box. 

The  electrochemical  measurements  were  performed  at  inser¬ 
tion  and  extraction  current  density  of  50  mA  g  1  between  voltage 
limits  of  0.01  V  and  ~3  V.  The  silicon  oxycarbide  materials  show  a 
first  cycle  loss  of  ~30%  but  are  stable  and  reversible  in  subsequent 
cycling.  Their  reversible  capacity  is  700—1000  mAh  g-1  at  about 
0.2  C  h  These  materials  can  sustain  very  high  C-rates  without 
failure  [8], 

3.2.  Galvanostatic  titrations 


5|r-T»  +  35e_20r(,-r2“”)  <20) 

Eq.  (20)  has  the  correct  long  time  limit.  That  is,  when  t*>t0,  the 
second  term  on  the  right  hand  side  vanishes,  giving  the  expected 
result  that  the  surface  concentration  becomes  equal  to  the  average 
concentration  within  the  particle  as  given  by  Eq.  (7). 

In  order  to  compare  Eq.  (20)  to  experiment,  it  is  useful  to  write  it 
in  the  form 


30t0 

c(t0)(1  -e- 

where 


-A ca{x*)-e  20t* 


(21) 


Aca  =  Caff)  -  C(T„). 

Eq.  (21 )  specifies  that  the  evolution  of  the  concentration  from  its 
initial  value,  ca(r*),  to  its  uniform  steady  state  value,  ca(ro),  is 
approximated  by  an  exponential  decay,  meaning  that  the  concen¬ 
tration  falls  to  (1/e)  of  its  initial  value  after  a  period  of  20t*. 


The  results  reported  here  are  in  the  stable  and  reversible  cycling 
regime.  The  starting  state  for  the  titration  experiments  consisted  of 
extracting  lithium  from  the  anode  with  a  current  of  50  mAh  g-1 
until  reaching  3  V,  at  which  point  the  anode  was  assumed  to  be 
essentially  depleted  of  Li. 

Current  was  then  injected  into  the  anode  (at  50  mAh  g-1)  for 
30  min,  equivalent  to  Li  insertion  of  25  mA  g_1.  Next,  the  current 
was  interrupted  and  the  voltage  was  measured  in  open  circuit, 
which  gradually  rose  with  time  towards  a  constant  value  (diffusion 
relaxes  the  surface  concentration  which  leads  to  a  higher  voltage). 
The  voltage  reached  a  steady  state  at  the  end  of  24  h;  at  this  point 
the  Li  concentration  in  the  particle  was  assumed  to  have  become 
uniform.  These  two  steps,  current  injection  followed  by  open  circuit 
relaxation,  constituted  one  episode  in  a  series  of  twenty  similar 
episodes  repeated  successively,  with  each  episode  building  on  the 
one  before.  The  results  are  given  in  Fig.  4.  The  data  are  plotted  in 
two  ways:  the  upper  graph  gives  the  cell  voltage  versus  time,  and 
the  lower  versus  capacity.  In  the  time  domain  the  current  injection 
portion  is  short  and  the  relaxation  portion  for  each  episode  is  long. 
In  the  capacity  domain  the  vertical  segments  represent  open  circuit 
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Step-wise  Insertion  of  Li  followed  by  Open-Circuit-Relaxation 


plotted  against  time  (upper  figure),  and  against  the  total  inserted  charge  (lower  figure). 
Adopted  from  Ref.  [1). 


relaxations  (at  constant  capacity)  which  are  followed  by  injection 
of  additional  Li  into  the  material. 

The  steady  state  values  of  the  voltage  from  each  of  the  open 
circuit  experiments  are  summarized  in  Table  1.  They  list  this  voltage 
and  the  corresponding  capacity  for  each  episode.  The  potential 


Steady  state  Nernst  potential  measured  after  open-circuit-relaxation  at  different 
levels  of  the  state-of-charge. 

Episode#  OCR  start  OCR  end  Capacity  Qu  Steady  state  State-of-charge 
potential 

h  h  mAh  g-1  V % 


changes  form  1.31  V  at  25  mA  g-1  to  0.16  V  at  652  mA  g-1.  The 
capacities  are  also  expressed  as  the  state-of-charge  (SOC)  by 
recognizing  that  the  capacity  when  charged  up  to  0.01  V  is 
1015  mA  g1.  The  start  and  end  time  for  the  open  circuit  relaxation 
for  each  episode  is  given  in  the  left  columns  of  Table  1. 

The  time  dependent  relaxations  contain  information  about  the 
diffusion  of  Li  from  the  surface  towards  the  center  of  the  particles. 
In  the  following  section  these  data  are  compared  with  theory  to 
obtain  a  value  of  the  chemical  diffusivity  of  Li  in  silicon-oxycarbide. 


4.  Determination  of  chemical  diffusivity 

We  seek  to  apply  the  theoretical  analysis  in  Section  2  to  the 
titration  results  reported  in  Section  3  in  order  to  elicit  a  value  for 
the  coefficient  of  chemical  diffusion.  This  is  a  two-step  procedure 
since  the  analysis  is  in  terms  of  the  concentration,  while  diffusion  is 
driven  by  gradients  of  the  chemical  potential.  The  relationship 
between  the  chemical  potential  and  the  concentration,  the  first 
step,  is  addressed  below. 


4.1.  Relating  chemical  potential  to  concentration 
The  transport  equation  for  chemical  diffusion  is 

Ju  =  (22) 

Here  D  is  the  coefficient  of  chemical  diffusion,  cu  is  the  volu¬ 
metric  concentration  of  the  diffusion  species  (Li),  with  units  of 
number  per  unit  volume,  whose  chemical  potential  is  written  in 
terms  of  its  activity,  an  as 

hu  =  p0u+RT\n(au).  (23) 

Here  the  standard  state  is  pure  Li  metal  with  an  activity  of  unity 
and  a  chemical  potential  of 

The  theoretical  analysis  in  Section  2  has  been  solved  for  Fickian 
diffusion  where  the  flux  is  proportional  to  the  gradient  of  the 
concentration  (in  units  of  number  of  species  per  unit  volume),  so 
that 

Ji  i  =  -DFVLi  (24) 

where  DF  is  the  Fickian  diffusion  coefficient  of  lithium  in  the  elec¬ 
trode  material. 

Taking  the  gradient  of  Eq.  (23),  substituting  into  Eq.  (22),  and 
then  equating  the  result  to  Eq.  (24)  gives  the  desired  relation  be¬ 
tween  D  and  Dp,  namely 


12 

13 

15 

16 
17 


20 


0.5  24.5  25 

25.0  37.0  50 

37.5  61.5  75 

62.0  86.0  100 

86.5  110.5  125 

111.0  135.0  150 

135.5  159.5  176 

184.5  208.5  201 

209.0  235.0  251 

235.5  257.5  276 

258.0  282.0  301 

282.5  306.5  326 

307.0  319.0  351 

320.0  332.0  401 

333.0  345.0  451 

346.0  358.0  502 

359.0  371.0  552 

372.0  384.0  602 

385.0  397.0  652 


1.31  2.2 

1.09  4.4 

0.96  6.4 

0.86  8.6 

0.78  10.8 

0.69  13.0 

0.66  15.3 

0.61  17.7 

0.52  22.6 

0.48  25.3 

0.45  28.0 

0.41  30.8 

0.37  34.8 

0.33  39.5 

0.28  44.5 

0.24  49.9 

0.21  55.0 

0.18  59.8 

0.16  64.2 


Dp  =  aD  (25) 

where 


=  dln(qLi) 
dln(cu) ' 

Furthermore,  combining  Eq.  (26)  with  Eq.  (2)  leads  to 


(26) 


and 


dV 

dlncLj ' 


(27) 


The  results  of  the  titration  experiments  given  in  Table  1  provide 
data  for  the  voltage  as  a  function  of  concentration  which  can  be 
used  to  obtain  the  value  for  ay. 

It  is  of  note  that  Eq.  (27)  may  be  used  to  calculate  the  concen¬ 
tration  from  the  measurement  of  the  voltage 
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£u 

4 


eV/ar 


(28) 


Table  2 

Analysis  of  the  open  circuit  relaxation  data  for  one  episode  to  obtain  the  relaxation 
time  which  is  then  matched  to  experiment  as  illustrated  in  Fig.  7. 


where  Cy  is  the  maximum  charge  held  by  the  anode  when  V  — »  0. 

One  still  needs  a  relationship  between  capacity,  measured  in 
units  of  mA  g_1,  and  Cu.  This  relation  becomes  redundant  because 
Eq.  (27)  is  written  in  terms  of  the  derivative  of  a  logarithm,  which  is 
a  dimensionless  quantity  evident  from 

dlncLi=^.  (29) 

cu 

We  denote  the  capacity  as  Qu  with  units  mA  g_1.  Since  the 
concentration  is  also  the  capacity  in  different  units,  Cu  can  be 
replaced  by  Ql,  to  obtain 


av 


dV 

(®S  Qu 


(30) 


Thus  the  experimental  value  for  av  can  be  determined  from  the 
slope  of  a  plot  of  V  versus  In  Qu.  Such  a  plot  of  the  data  given  in 
Table  1  is  shown  in  Fig.  5.  The  data  fit  a  straight  line  whose  slope 
gives  ay  =  0.35  V.  Substituting  this  value  into  Eq.  (27)  gives  a  —  13.5. 

A  constant  value  for  a,  arising  from  a  straight  line  fit  to  the  data 
through  the  full  range  of  capacity,  gives  a  linear  proportionality 
between  the  Fickian  and  the  chemical  diffusion  coefficients,  which 
greatly  simplifies  the  application  of  the  analysis  to  the  experiment. 
The  linear  fit  is  ascribed  to  the  amorphous  nature  of  the  silicon- 
oxycarbides  and  the  absence  of  phase  change  with  Li  capacity. 
The  case  of  systems  that  encounter  phase  change  is  discussed  in 
Section  6. 

The  next  step  is  to  match  the  time  dependent  relaxation  of  the 
voltage  under  open-circuit  to  Eq.  (21 )  in  order  to  extract  a  value  of 
the  time  constant  At*,  in  real  time  which  then  gives  the  value  for  DF. 


4.2.  Diffusion  during  open  circuit  relaxations 

During  open  circuit  relaxation  Li  diffuses  within  the  particle 
toward  its  terminal  uniform  concentration.  Theoretically,  the 
relaxation  is  given  by  Eq.  (21 ).  The  left  hand  side  in  this  equation  is 
the  difference  between  the  surface  concentration  at  time  t*  and  the 
steady  state  concentration  achieved  after  full  relaxation.  The 
additional  factor  30to/c(to)(1  -  e~20l°)  on  the  left  is  a  constant  that 
depends  on  the  magnitude  and  duration  of  current  injection,  io;  it 
is  independent  of  time  and  therefore  does  not  affect  the  mea¬ 
surement  of  the  time  constant  for  relaxation. 

The  problem  is  now  well  posed  since  the  time  dependent  con¬ 
centration  can  be  deduced  from  the  measurement  of  the  voltage 


ln{QLi ,  mAhg-1) 


Fig.  5.  Plot  of  data  according  to  Eq.  (27)  for  the  determination  of  av. 


Relative  time  Measured  Concentration  Relative  Percent  relative 

for  OCR,  f  cell  voltage  (arb.  u.)  concentration  concentration 


0.00  -1.2912 

0.02  -1.2963 

0.12  -1.3075 

0.22  -1.3134 

0.32  -1.3169 

0.46  -1.3202 

0.62  -1.3225 

0.86  -1.3248 

1.12  -1.3263 

1.39  -1.3278 

1.69  -1.3286 

1.99  -1.3291 

2.39  -1.3301 

2.82  -1.3306 

3.19  -1.3311 

3.56  -1.3312 

3.89  -1.3314 

4.29  -1.3316 

4.69  -1.3317 

5.09  -1.3317 

Steady  state  after  24  h 
24.0  -1.3358 


0.024993  0.002993 

0.024632  0.002632 

0.023856  0.001856 

0.023457  0.001457 

0.023224  0.001224 

0.023006  0.001006 

0.022855  0.000855 

0.022706  0.000706 

0.022609  0.000609 

0.022512  0.000512 

0.022460  0.000460 

0.022428  0.000428 

0.022364  0.000364 

0.022332  0.000332 

0.022301  0.000301 

0.022294  0.000294 

0.022281  0.000281 

0.022269  0.000269 

0.022262  0.000262 

0.022262  0.000262 


0.022003  0.000003 


40.9 

33.6 

28.6 
23.6 

20.3 

17.1 

15.4 
14.3 

12.2 


10.0 

9.8 
9.4 
9.0 

8.8 
8.8 


100 


during  the  relaxation.  The  relaxation  data  were  matched  to  Eq.  (21 ) 
for  each  of  the  twenty  episodes  at  different  values  of  the  state-of- 
charge.  The  method  is  outlined  in  Table  2  using  one  episode  as  an 
example.  The  first  column  contains  the  (real)  time  from  the  start  of 
the  open-circuit  experiment.  The  second  column  gives  the  value  of 
the  measured  cell  voltage.  The  third  column  converts  the  voltage  to 
the  dimensionless  concentration  relative  to  the  steady  state  value 
obtained  after  24  h,  which  is  given  at  the  bottom  of  the  table.  The 
last  column  reduces  the  concentration  to  a  fractional  value  starting 
with  100%  at  time  zero  and  approaching  zero  after  24  h. 

The  next  step  is  to  match  the  time  dependent  decay  in  con¬ 
centration  given  in  Table  2  to  Eq.  (21).  The  procedure  is  shown  in 
Fig.  6.  The  main  graph  shows  the  time  dependent  decay  in  con¬ 
centration  from  the  right  hand  columns  in  Table  2.  The  objective  is 


20t*  =  1 

Fig.  6.  Procedure  for  determination  of  the  relaxation  time  constant  by  fitting  a  single 
exponential  term  to  the  data  as  shown  in  the  inset. 
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to  match  this  graph  to  Eq.  (21 ).  The  left  hand  side  of  Eq.  (21 )  is 
dimensionless.  It  represents  the  change  in  concentration  from  100% 
at  r*  =0  to  zero  as  i*  ->  °° .  We  give  it  the  simpler  form 

y  =  e“20T*  (31) 

where  y  represents  the  dimensionless  value  of  the  surface  con¬ 
centration.  Eq.  (31 )  is  plotted  in  the  inset  of  Fig.  6.  When  overlaid  on 
the  main  graph  it  provides  the  equivalence  between  the  dimen¬ 
sionless  and  real  time.  The  equivalence  between  the  dimensionless 
time  constant  20r*,  and  the  real  time  constant  At*,  is  now  self 
evident.  Thus  when  20r*  =  1,  the  real  time  is  equal  to  At*  =  0.3  h. 

The  definition  of  normalized  time  given  in  Eq.  (5)  leads  to  the 
diffusion  coefficient  by  rewriting  it  in  the  form 

20r*  =  20  =  1  (32) 

az 

giving 


10-15  to  10-10  m2  s_1  [9],  ft  is  not  clear  if  this  difference  is  actual,  or 
whether  it  can  be  attributed  to  the  different  ways  in  which  the 
diffusion  coefficients  have  been  measured.  The  wide  variability  in 
the  data  for  graphite  in  the  literature  may  reflect  uncertainties  in 
the  measurement  methods.  The  narrow  band  for  the  diffusion 
coefficient  measured  over  a  wide  range  of  the  state-of-charge  is  a 
notable  feature  of  the  present  work. 

4.3.  The  in  fluence  of  particle  size 

In  Eq.  (33)  the  diffusion  coefficient  is  a  material  parameter. 
Therefore  the  time  constant  for  relaxation,  At*,  is  related  to  the 
particle  radius,  a.  We  note  in  Fig.  7  that  the  time  constant  for 
relaxation  is  quicker,  about  0.3  h  for  SOC  =  2.2%.  This  is  readily 
explained  by  the  early  saturation  of  small  particles  in  the  powder.  It 
can  be  inferred  that  the  2.2%  of  the  SOC  is  held  by  these  particles.  It 
follows  that  the  volume  fraction  of  these  small  particles  in  the 
powder  was  also  about  2.2%.  A  scanning  electron  micrograph  of  the 
silicon-oxycarbide  powder  is  given  in  Fig.  8.  Note  that  while  most 
particles  have  a  radius  of  approximately  5  pm,  there  is  a  sprinkling 
of  particles  that  are  much  smaller,  which  we  hold  responsible  for 
the  short  relaxation  time. 


where  for  the  episode  illustrated  in  Fig.  6,  At*  =  0.3  h.  Substituting 
this  value  and  the  particle  size  into  Eq.  (33)  gives  the  value  for  D f. 

Each  of  the  twenty  relaxations  were  analyzed  using  the  proce¬ 
dure  outlined  above.  They  cover  a  wide  range  of  the  state-of-charge 
(SOC),  from  2.2%  to  64.2%.  These  results  are  shown  together  in  Fig.  7. 
Remarkably,  all  relaxations  ranging  form  4.4%-64.2%  SOC  fall 
within  a  relatively  narrow  band  of  At*  =  1.5  -  3.5  h.  Considering 
that  the  SOC  changes  by  more  than  a  factor  of  10,  this  result  reflects 
good  support  for  the  main  result  of  the  analysis:  the  relaxation  can 
be  expected  to  be  independent  of  the  duration  and  extent  of  the 
charge  injection  into  the  electrode  particles. 

For  the  case  of  SOC  =  4.4%-64.2%,  we  assume  a  particle  radius 
a  ~  2.5  pm.  Substituting  At*  =  1.5  -  3.5  h  into  Eq.  (33)  gives 
Dp  =  2.5  x  10-17— 5.8  x  10-17  m2  s-1.  Further  substitution  into  Eq. 
(25),  a  =  13.5,  gives  the  following  value  for  the  chemical  diffusion 
coefficient  of  Li  in  silicon-oxycarbide 

D  =  1.8  x  10  18— 4.2  x  10~18  m2  s_1  (34) 

The  value  in  Eq.  (34)  is  lower  than  the  values  for  diffusion  of  Li  in 
graphitic  structures  reported  in  the  literature,  which  range  from 


t  *,  h 


Fig.  7.  Plots  of  surface  concentrations  in  real  time  for  all  20  episodes  of  open-circuit 
relaxations.  Each  plot  followed  the  procedure  outlined  in  Table  2. 


4.4.  Analysis  of  current  injection  data 

The  plots  in  Fig.  3  describe  the  results  for  the  change  in  surface 
concentration  as  the  current  is  injected  into  spherical  particles  at  a 
constant  rate.  The  injection  time  is  prescribed  by  t,  a  normalized 
parameter,  which  is  related  to  real  time  via 


The  above  equation  is  obtained  by  combining  Eq.  (5)  with  Eq. 
(33).  Here  t  is  real  time,  and  At*  is  measured  from  the  open-circuit- 
relaxation  experiments. 

It  now  becomes  possible  to  compare  the  change  in  the  surface 
concentration  during  current  injection,  as  predicted  in  Fig.  3,  with 
experiment.  This  comparison,  for  episode  one,  is  shown  in  Fig.  9. 
Here  the  time  axes  for  theory  and  experiment  have  been  syn¬ 
chronized  by  setting  At*  =  0.3  h,  the  value  measured  for  episode 
one.  For  example,  r  =  0.08  in  dimensionless  time  is  equivalent  to 
0.5  h  in  real  time.  The  vertical  axis  has  arbitrary  units.  The  result  in 
Fig.  9  suggests  that  the  solid-electrolyte-interface  and  the  ionic 
resistance  of  the  electrolyte  do  not  significantly  affect  the  measured 
cell  voltage  for  the  silicon-oxycarbide  materials  when  the  current 
density  is  50  mAh  g  This  may  not  necessarily  be  the  case  at 
higher  current  densities. 


5(im 

Fig.  8.  Scanning  electron  micrograph  of  silicon-oxycarbide  particles. 
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0  0.25  0.5 


T 


Fig.  9.  Comparison  of  experiment  and  theory  for  the  change  in  surface  concentration 
during  current  injection.  The  normalized  time  and  real  time  have  been  synchronized 
along  the  abscissa. 


5.  Discussion 

In  this  section  we  compare  the  present  results  to  similar  work  in 
the  literature.  They  are  summarized  below.  The  current  analysis 
and  experiments  are  concerned  with  galvanostatic  titrations  where 
intermittent  injection  of  charge  at  constant  rate  is  followed  by 
measurement  of  the  time  dependent  relaxation  of  the  galvanic 
potential  at  open  circuit. 

(i)  Effect  of  particle  size 

According  to  Eq.  (33)  the  diffusion  coefficient  is  related  to  both 
the  relaxation  time  and  the  particle  radius.  While  the  scaling  in  this 
result  is  consistent  with  Eq.  ( 1 ),  it  differs  from  it  by  a  large  nu¬ 
merical  factor.  The  measurements  of  the  relaxation  time,  At*,  as  a 
function  of  the  state-of-charge  given  in  Fig.  7  highlight  how  its 
value  can  vary  because  small  particles  become  “saturated”  with  Li 
quickly.  (The  small  particles  normally  occupy  minor  volume  frac¬ 
tions  since  volume  is  proportional  to  the  third  power  of  particle 
size.)  This  result  is  significant  for  understanding  how  the  particle 
size  distribution  can  affect  the  C-rate  capability  of  electrode 
materials. 

(ii)  The  three-dimensional  nature  of  the  problem 

The  conventional  method  for  deriving  the  diffusion  coefficient 
from  time  dependent  relaxations  has  been  to  consider  one¬ 
dimensional  diffusion  in  a  semi-infinite  slab  [3-5],  which  gives 
solutions  in  the  form  of  error  functions  [10],  Here  an  exact  solu¬ 
tion  for  the  three-dimensional  problem  is  presented,  and  it  is 
rigorously  shown  that  the  relaxation  time,  to  a  very  good 
approximation,  can  be  described  by  a  single  exponential,  given  by 
Eq.  (31). 

The  spherical  geometry  has  also  been  considered  in  recent  work 
[6,11  ]  where  the  time  dependent  problem  is  approached,  not  in  the 
real-time,  but  in  the  Laplace  transform  domain.  In  this  approach  it 
is  more  difficult  to  extract  non-dimensional  parameters  which  can 
be  unambiguously  and  critically  related  to  experiments.  The  com¬ 
parison  between  theory  and  experiment  in  Figs.  6,  7  and  9,  are 
examples  of  the  real-time  approach:  in  these  graphs  only  the  time 
constants  for  relaxation  are  needed  since  the  concentrations  can  be 
allowed  to  have  arbitrary  units. 


(iii)  Achieving  convergence  between  theory  and  experiment 

We  have  assiduously  sought  analytical  results  that  can  be 
applied  to  experiments,  by  employing  non-dimensional  parame¬ 
ters.  For  example  it  is  noteworthy  that  the  relaxation,  given  by  the 
right  hand  side  of  Eq.  (21 ),  is  independent  of  the  state-of-charge, 
and  also  independent  of  the  rate  of  insertion  of  the  charge.  The 
only  condition  is  that  the  relaxation  be  measured  from  the  moment 
that  the  current  stops  and  the  open-circuit  condition  begins.  This  is 
a  significant  simplification  in  the  application  of  the  analysis  to 
experiment.  Thus,  the  usual  caveat  that  the  current  must  be 
injected  very  quickly  in  the  GITT  method  is  not  necessary. 

(iv)  Experimental  determination  of  chemical  activity  as  a  func¬ 
tion  of  concentration 

In  order  to  compare  theory  with  experiment,  it  is  necessary  to 
convert  the  cell  voltage  into  Li  concentration,  which  requires 
knowledge  of  ay,  which  is  defined  in  Eq.  (30).  Its  value  is  obtained 
by  plotting  the  logarithm  of  the  state-of-charge  against  the  satu¬ 
ration  value  of  Nernst  voltage  measured  in  open  circuit.  Caution 
must  be  exercised  if  ay  can  change  should  the  material 
undergo  phase  transition.  This  point  is  discussed  further  in  the  next 
section. 

6.  The  nuances  of  theory  and  experiment 

Engineering  materials,  and  the  systems  that  employ  them,  are 
nearly  always  complex.  For  example,  the  performance  of  anodes  in 
a  cell  depends  on  the  morphology  and  microstructure  of  the  active 
material,  the  architecture  of  the  electrode,  and  the  various  kinetic 
and  interfacial  reactions. 

Efforts  to  model  such  engineering  systems  often  follow  two 
approaches.  One  is  to  employ  numerical,  computational  methods 
that  seek  to  simulate  the  full  complexity  of  the  system.  The 
objective  of  these  simulations  is  to  calculate  the  numerical  values  of 
the  experimental  outcomes.  These  simulations  often  do  not  provide 
simple  parametric  relationships  between  material  parameters  and 
cell  performance. 

In  the  second  approach  the  phenomenological  information  from 
experiments  is  distilled  into  concepts  that  are  developed  into 
analytical  models.  The  objective  of  these  closed-form  models  is  to 
relate  the  microstructural  parameters  to  experiment.  The  models 
build  upon  fundamental  relationships  from  basic  science.  While 
these  comparisons  require  approximations,  they  serve  to  provide 
scaling  relationships  between  microstructure  and  experimental 
variables.  The  deviation  of  experiment  from  theory  can  provide 
new  insights  into  atomistic  mechanisms.  In  other  words  this 
approach  enables  transparency  between  science  and  engineering. 

The  closed  form  models  necessarily  require  idealized  geome¬ 
tries.  For  example,  the  actual  shape  of  the  particles,  shown  in  Fig.  8 
is  clearly  not  spherical  as  has  been  assumed.  These  assumptions 
require  a  nuanced  approach  to  comparing  theory  with  experiment, 
which  we  discuss  below. 

We  address  three  issues:  (i)  the  cuboid  shape  of  the  particles  as 
in  Fig.  8,  (ii)  the  influence  of  phase  transitions  when  Li  is  inserted 
into  an  anode  material,  and  (iii)  the  separation  of  the  time  con¬ 
stants  during  relaxation  to  determine  which  one  is  most  likely  to  be 
related  to  the  diffusion  of  Li  into  the  volume  of  the  anode  material. 

The  error  from  the  shape  of  the  particles  is  addressed  by 
establishing  the  extent  to  which  the  shape  deviates  from  spherical. 
Exact  results  for  sphere  and  cuboid  shapes,  obtained  by  FEM  pro¬ 
vide  an  estimate  of  this  error.  An  optical  micrograph  of  the  particles 
taken  in  a  cross-sectional  view  is  given  in  Fig.  10.  The  particles 
are  angular  and  look  approximately  like  cuboids  with  an  aspect 
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Fig.  10.  Optical  micrograph  of  a  cross-sectional  view  of  silicone  oxycarbide  particles. 


ratio  (the  smallest  to  the  largest  dimension)  in  the  range  of  about 
0.5-0.75. 

Exact  results  for  sphere  and  cuboid  shapes,  obtained  by  finite- 
element-analysis  (FEM)  provide  an  estimate  to  this  error.  To 
determine  the  extent  to  which  the  results  for  the  sphere  would  be 
representative  for  these  shapes,  we  have  carried  out  numerical  FEM 
analysis  to  determine  the  time  for  surface  concentration  to  reach 
the  interior  of  the  particle.  The  results  for  a  cuboid  with  an  aspect 
ratio  of  0.5,  a  cube,  with  an  aspect  ratio  of  1.0,  and  a  sphere  all  with 
the  same  width  are  shown  in  Fig.  11.  In  each  case  a  concentration  of 
unity  was  applied  to  the  surface  as  a  step  function,  and  the  time- 
dependent  concentration  at  the  center  of  the  body  was  calcu¬ 
lated.  The  cube  is  approximately  two  times  slower  to  reach  the 
same  concentration  at  its  center  than  the  sphere.  The  sphere  and 
the  cuboid  are  very  similar.  The  micrograph  in  Fig.  10  shows  that 
most  of  the  particles  have  an  aspect  ratio  of  approximately  0.5.  The 
inference  is  that  the  sphere  was  quite  a  good  approximation  for  the 
diffusion  analysis.  Even  assuming  the  shapes  to  be  cubes  would 
lead  to  an  error  of  no  more  than  a  factor  of  two  for  the  diffusion 
coefficient. 

The  coefficient  of  chemical  diffusion,  D,  is  proportional  to  the 
Fickian  diffusion  coefficient,  Df,  through  a,  as  given  in  Eq.  (25). 
Experiments  with  silicon-oxycarbide  (SiCO)  have  shown  that  a  has 
a  constant  value.  This  is  perhaps  a  special  case  since  SiCO  is 


sphere 


Relative  Time 

Fig.  11.  Plots  of  the  concentration  at  the  core  of  three  particle  shapes  with  time,  after 
the  step-function  application  of  unit  concentrations  at  the  surface  at  time  zero.  The 
results  were  obtained  from  numerical  finite  element  analysis.  Note  that  the  cuboid 
with  an  aspect  ratio  of  0.5  behaves  similarly  to  the  sphere.  The  perfect  cube  is  slower 
by  a  factor  of  about  two. 


In  (lithium  concentration,  a.u.*) 


Fig.  12.  Conceptual  plots  of  Eq.  (A7)  in  a  crystalline  material  A  which  transforms  to  a 
new  phase  B  when  Li  is  inserted  into  it.  Distinct  slopes  for  uv  can  be  determined  only  in 
the  single  phase  regimes. 


amorphous  and  remains  so  over  the  full  range  of  the  state  of  charge. 
Crystalline  materials,  however,  often  encounter  a  phase  transition 
with  Li  concentration.  Let  us  say  these  phases  are  A  and  B,  with  the 
first  being  stable  at  low  and  the  second  at  high  Li  concentrations. 
However,  there  would  also  be  a  range  where  both  phases  exist  in 
equilibrium  with  Li  partitioned  between  them  as  required  by 
thermodynamic  equilibrium.  In  this  two-phase  region,  Li  is  ex¬ 
pected  to  partition  increasingly  to  the  second  phase  at  the  expense 
of  the  first,  at  constant  voltage.  In  this  two-phase  region  the 
diffusion  of  Li  will  be  complex.  Experimental  plots  of  the  type 
shown  in  Fig.  12  may  be  used  to  identify  the  single  phase  regions, 
and  the  values  for  and  a®  for  each  phase  determined  from  the 
slopes  in  the  single  phase  regions. 

The  third  question  pertains  to  the  irreversible  capacity  often 
seen  in  the  first  cycle.  The  experiments  and  the  theory  in  this  article 
refer  to  reversible  insertion  of  lithium,  with  nearly  full  coulombic 
efficiency.  The  irreversible  capacity  is  sequestered  into  deep  energy 
levels  where  it  can  be  assumed  to  be  immobile  under  the  nominal 
range  of  voltage  levels.  One  of  the  key  results  from  the  theoretical 
analysis,  which  is  borne  out  by  experiments,  is  that  the  relaxation 
profile  is  independent  of  the  prior  history  of  current  injection. 
Furthermore,  the  relaxation  time  is  independent  of  the  rate  and  the 
time  period  of  current  injection.  The  only  condition  is  that  the  Li 
concentration  in  the  particle  be  uniform  before  current  injection. 
This  feature  of  the  results  lends  considerable  generality  to  the 
galvanostatic  titrations. 

Finally,  we  consider  the  possible  role  of  extraneous  relaxations  on 
the  measurement  of  the  diffusion.  Equilibration  of  Li+  within  the 
electrolyte,  and  charge  transfer  at  anode-electrolyte  interfaces  are 
two  examples  of  such  effects.  Fortunately,  each  effect  has  its  own  time 
constant  for  relaxation.  We  do  see  very  short  relaxations  at  the  onset 
of  the  open  circuit  relaxation.  They  have  negligibly  small  voltage 
amplitudes.  The  time  constants  for  these  extraneous  relaxations  are 
typically  less  than  1  min.  By  comparison,  the  Li  diffusion  time  con¬ 
stants  approach  an  hour.  In  our  experiments  the  effect  of  the  extra¬ 
neous  relaxations  has  been  filtered  out  by  taking  data  5  min  after 
initiating  the  open  circuit.  Further  support  for  the  diffusive  relaxations 
lie  in  their  exponential  behavior,  implying  a  single  mechanism.  The 
systematic  and  large  voltage  amplitude  of  the  relaxations  further 
underpins  the  interpretation  of  Li  diffusion  within  the  particles. 

7.  Summary 

1.  Open  circuit  relaxations,  following  stepwise  injections  of  Li 

currents  into  spherical  particles  are  analyzed.  The  relaxations 
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can  be  described  by  Eq.  (21)  which  reduces  to  Eq.  (31)  for 
application  to  data.  Thus,  the  measurements  for  different  levels 
of  injection  rate,  and  state-of-charge,  can  be  represented  in  a 
universal  curve.  That  the  relaxation  time  constants  are  inde¬ 
pendent  of  the  rate  of  current  injection  and  the  state  of  charge 
imparts  great  generality  and  simplicity  to  the  galvanostatic 
titration  experiments.  In  the  past  it  has  been  assumed  that  for 
these  experiments  to  be  valid  the  current  must  be  injected  in  a 
short  pulse.  We  show  that  this  condition  is  unnecessary. 

2.  Steady  state  Nernst  potentials  measured  at  different  states  of 
charge  are  plotted  according  to  Eq.  (27)  to  obtain  a  value  for  ay. 
Then,  Eq.  (28)  can  be  used  to  convert  the  voltage  into  concen¬ 
tration  at  particle  surfaces  during  open  circuit  relaxations. 

3.  The  relaxation  time  At*  is  related  to  the  Fickian  diffusion  coef¬ 
ficient  through  Eq.  (33). 

4.  A  distribution  in  the  particle  size  can  produce  a  spread  in  the 
relaxation  time.  The  small  particles,  which  would  have  a  shorter 
relaxation  time,  can  saturate  at  low  values  of  the  state-of-charge. 
Therefore  the  relaxation  time  can  increase  with  increasing  state- 
of-charge  if  the  particles  have  a  broad  size  distribution. 

5.  The  change  in  the  surface  concentration  of  the  spherical  parti¬ 
cles  during  charge  injection  is  analyzed.  This  result  is  given  in 
Eq.  (11 ).  A  comparison  of  experiment  with  this  analysis  provides 
insight  into  the  significance  of  interfacial  polarization  and 
resistive  potential  drops  across  the  electrolyte. 
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Appendix  A.  Constant  current  solution 

1  Introduction 


^  J  c(r,t)r2dr  =  a2J0. 
o 

Integrating  on  time  we  find 


J  c(r,t)r2 dr  =  a2J0t. 
o 

The  volume-averaged  concentration  is  defined  as 


1=1  |  c(r,t)dV  =  A  J  C(r,t)r2dr 


so  that  the  volume-averaged  concentration  grows  linearly 


c  =  a^  =  3J°t. 
a3/  3  a 


(A4) 


(A5) 


(A6) 


(A7) 


1.2  Transient  solution 

Inserting  (A7)  into  (Al)  shows  that  the  transient  solution  will 
have  a  term  depending  on  r,  thus  we  write 

c(r,t)  =  c(t)+F(r)  +  v(r,t).  (A8) 

into  the  diffusion  equation  gives 


and  the  boundary- value  problem  for  v(r,t)  becomes 


(A9) 


(AlOa) 


In  this  regime  for  which  0  <  t  <  to  the  current  flows  from  the 
anode  to  the  cathode  resulting  in  a  flux  Jo  of  lithium  ions  into  the  v(r,0)  =  -^,  =0, 

sphere.  Here  one  must  solve  the  axisymmetric  diffusion  equation  in  2aDp  dr  |r=0  dr  \r=a 

spherical  coordinates,  viz. 


0. 

(AlOb) 


dc 

at 


2  3c\ 
FarJ- 


(Al) 


The  constant  flux  of  ions  through  the  surface  of  the  sphere  and 
the  initial  concentration  of  ions  being  zero  lead  to  the  initial  and 
boundary  conditions 


TL,^'  <«> 

Since  flux  is  continuously  added  over  time,  the  sphere  contin¬ 
uously  increases  in  ion  concentration  without  bound.  This  is  shown 
in  the  following  section. 


1.1  Volume  averaged  concentration 

We  integrate  conservation  Eq.  (Al )  over  the  volume  to  obtain 


(A3) 


Now  the  problem  is  an  initial  distribution  of  concentration 
with  zero  flux  through  the  sphere.  To  solve  this  we  transform  ac¬ 
cording  to 

w(r,  t)  =  rv(r,  t)  (Al  1 ) 

to  obtain 


aw  _  a2w 

aF  “  V 


(A12a) 


WP-O, 

(A12b) 

The  Robin  boundary  condition  causes  no  problem  because  it  is 
homogeneous  and  we  therefore  have  homogeneous  boundary 
conditions  with  an  initial  distribution  of  ‘concentration’  in  these 
variables.  Assuming  the  separable  solution 


which  gives 


w(r,  t)  =  F{r)G{t) 


(A13) 
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we  find 

F"  +  ^F  =  0,  F(  0)  =  0,  F'(a)  =  ^F(a)  (A14) 

C'  +  AG=0,  G(0)  =  1  (A15) 

where  A  is  the  separation  constant  and  primes  denote  differentia¬ 
tion  with  respect  to  the  arguments  of  F  and  G.  The  solutions  satis¬ 
fying  F(0)  =  0  and  G(0)  =  1  are 

F(r)  =  sin  a„r,  G(t)  =  e^D^t  (A16) 

where  the  separation  constant  has  been  set  to  A  =  DF a2.  Satisfying 
the  Robin  boundary  condition  in  (A14)  gives  the  eigenvalue 
equation 

/?„  =  tan/J„  (A17) 

where  f3n  =  aa„.  Inspection  of  Eq.  (A17)  shows  that  there  is  an 
additional  eigenfunction  at  A  =  0  given  by  F[r)  =  r  so  we  summarize 
the  eigenfunctions  as  follows 

</>o(r)  =  r,  (n  =  0)  (A18a) 

0n(r)  =  sinanr,  (n>l).  (A18b) 

Application  of  standard  Fourier  analysis  gives  the  solution 


w(r,  t)  =  (A19) 

o 

To  satisfy  the  initial  condition  we  must  have 


w(r,  0)  =  E  C„0n(r)  =*  (A20) 

for  which  the  constants  Co  and  C„  for  (n  >  1)  are  readily  found  to  be 


=  3  J0a  2J0 

0  10  DF’  "  DFa2sin(aan) 


(A21) 


Use  of  eigenvalue  Eq.  (A18)  has  been  made  to  obtain  the 
simplified  form  for  Cn  given  above.  Now  the  solution  for  w(r,t)  is 


w(r, t}  =  ^10Dp  r  +  £  Cnsin(anr)e  (A22) 

and  thus 

y(r, t )  =  +  1  £  CnsinCanrJe-^f.  (A23) 

Finally,  from  Eq.  (A8)  we  obtain,  on  reversion  to  fin,  the  final 
solution 

M  3J0  „  ,  J0(5r2  -  3a2)  2 J0a2  1  ^  sin(^nr/a)e-D^t/fl2 

_  a  t+  lOaDjj  DF  rZ. 

(A24) 

In  terms  of  dimensionless  radius  R  =  r/a  and  time  t  =  Dpt/a2  we 


C(*.T)  3t  1 

aJ0/DF  10 


(5s2-3)-i? 


sin 

^ns*n  fin 


(A25) 


The  normalized  surface  values  can  be  easily  computed  as 
function  dimensionless  time  t  by  evaluating  (A25)  at  R  =  1  which 
gives 


_  3T  _  I _  2 ye  ^ 
q/o/Dp  5  ^  ^  • 


(A26) 


Appendix  B.  Open  circuit  relaxation 

I  Introduction 


In  this  case  the  ion  current  to  the  sphere  stops  at  t  =  to,  and  the 
concentration  gradient  within  the  sphere  relaxes  by  diffusion  to  a 
final  uniform  concentration.  Being  interested  now  in  times  greater 
than  to  we  introduce  the  new  time  variable  t*  =  t  -  to. 


1.1  Final  steady  state 

Before  determining  the  time-dependent  adjustment  to  the  final 
uniform  concentration  we  compute  the  final  steady-state,  volume- 
averaged  concentration.  To  this  end  one  simply  evaluates  (A24)  at 
t  =  to  and  then  computes  the  volume  averaged  concentration  as 

Ju(r,t0)dV  (Bl) 


In  performing  the  integration  it  is  found  that  contributions  from 
the  r-dependent  terms  in  (A24)  are  zero,  so  the  only  contribution  is 
from  the  first  term  which  gives  the  final  equilibrated  concentration  as 


u(t- co  )  ^ 


3Jq[q 

a 


(B2) 


1.2  Adjustment  to  final  steady  state 

We  now  solve  the  diffusion  equation  given  (Al )  with  no  flux  at 
the  sphere  wall  and  with  the  known  distribution  of  concentration 
calculated  from  (A24)  at  the  stopping  time  t  =  to  for  which  t*  =  0. 
Thus  the  boundary  conditions  for  (Al )  are 


c(a,t*)  =  0, 


dc(r,  0)1 
dr  H4 


c(r,  0)  =  fir)  =A  +  Br2  -  E 


(B3a) 


where 

a  _  3/o  <0  _  3/qQ  jy  _  Jq  _  _  2Jpa _  Dp^to/a2 

a  10Df’  —  2aDp’  m  ~  D^2msinfim 

(B3b) 

Again  we  transform  according  to 

v(r,t*)  =  rc(r,  t*)  (B4) 

to  obtain  the  initial  and  boundary-value  problem 

W-H'S-  ',i0’n-0’  “tfE  = 

(B5) 
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From  the  analysis  in  Appendix  A  we  know  that  the  solution  is 
given  by 


C  =  2  J  -  2-fagZ  1  (\  _  e-Dfftto/a2' 

"  a^sin/3„  df  ft,  ^ 


v(r,t*)  =  y"Cn0n(r)e  DF“"r.  f}n  =  tan  (3n.  0n  =  aan.  (B6a) 


Now  we  can  compute  the  subsidiary  solution 


where 

0o(r)  =  r,  0„(r)  =  sin(a„r)  (B6b) 

We  now  determine  the  coefficients  C„.  The  first  coefficient  Co  is 
given  by 

J  0o(r)r/(r)dr 

C0  =  # - - - =  (B7) 

J  0o(r) 

o 

and  we  now  must  calculate 


v(r,t*)  =  C0Mr)  +  £  C„0n(r)e-^^r 


and  apply  (B4)  to  obtain 

c(r,  f)  =  3Jot°  +  2“io  V  ^r1  (l  - 

a  D  n^\P2sin0„  '  ' 

sin  (Pnr/a)  thfnr/<P 

In  dimensionless  form  we  have 


ah/D  0  r2 r  ’^sin  fin  ‘ 


(BIO) 


(Bll) 


(B12) 


(B13) 


/“  where  z*  —  Dpt* /a2. 

rsin(a„r)f  (r)dr=J  (B8a)  One  can  also  determine  the  time  evolution  of  the  surface  con- 

o  centration  for  this  initial  value  by  setting  R  =  1  in  (B13)  to  obtain 


where 

J=h  +J2-J3 


aJo/Dp 


-  3t0  =4 


2E^- 


(B14) 


h=AJ  rsin(anr)dr  =0,  J2  =  B  J  r3sin(anr)dr 
o  o 

 2Ba4sin 


For  J3  we  have 

h  =  E  pk  J  sin(akr)sin(anr)dr  =  °P„sin2  0n.  (B9) 

(f=t  0 

Then  one  obtains  J  =  ]2  -  h  and  the  coefficients  for  n  >  1  are 
given  by 
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